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PREFACE 

The  determination  of  diffusely  transmitted  and  re¬ 
flected  radiation  fields  for  an  atmosphere  of  finite  thick¬ 
ness  is  a  basic  problem  in  radiative  tr-rsfer.  This  Memoran¬ 
dum  presents  methods  of  calculating  these  fields  for  a  model 
planetary  atmosphere,  where  the  surface  of  the  planet  is 
assumed  to  scatter  incident  radiation  according  to  Lambert's 
law  of  scattering.  Scattering  in  the  atmosphere  is  assumed 
to  be  isotropic,  and  curvature,  polarization,  and  frequency 
dependence  are  not  considered.  The  Memorandum  also  provides 
numerical  results  and  checks. 


SUMMARY 


This  Memorandum  provides  formulas  for  obtaining  the 
diffusely  transmitted  and  reflected  radiation  fields  for 
a  planetary,  isotropically  scatter in_  atmosphere  of  finite 
thickness  in  terms  of  the  solution  to  the  problem  with  no 
planetary  surface.  Numerical  results  show  that  these  re¬ 
flected  and  transmitted  fluxes  are  essentially  the  same 
whether  isotropic  or  Rayleigh  scatter’ rg  laws  are  assumed. 
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I.  INTRODUCTION 

A  fundamental  problem  in  radiative  transfer  studies 
of  a  finite  atmosphere  illuminated  by  sunlight  is  the 
determination  of  the  intensities  of  the  diffusely  reflected 
and  transmitted  radiation  (i.e.,  the  radiation  reaching  the 
surface).  In  a  recent  paper  [1],  Kahle  studied  this  problem 
for  the  case  of  a  plane-parallel  atmosphere  bounded  by  a 
reflecting  surface  obeying  a  Lambert  law  of  reflection.  A 
Rayleigh  law  of  scattering  was  assumed  in  the  atmosphere, 
and  reflection  and  transmission  functions  were  obtained  via 
the  or 1 '  *- ;  ~~  "f  fhe  appropr ;  ate  sir.gn~ur  integral  eauations. 

This  Memorandum  presents  an  alternate  anproach  to  a 
related  planetary  problem.  Isotropic  rather  than  Ravleigh 
scattering  is  assumed  in  the  atmosphere.  Aside  from  the 
intrinsic  interest  of  the  problem,  a  maier  goal  of  this 
work  is  to  ascertain  the  differences  in  reflected  and  trans¬ 
mitted  fluxes  caused  by  these  two  scatter i nu  laws.  This 
Memorandum's  conclusion  is  that  virtually  no  difference 
exists . 
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II.  DERI  VAT  101.  OF  THE  EQUATIONS 


To  derive  an  apnropriate  system  of  ordinary  differ¬ 
ential  equations  for  the  reflected  and  transmitted  inten¬ 
sities  in  the  case  of  a  plane-parallel  atmosphere,  define 
the  function 


r(v,u,x)  =  R(v,u,x)/4v 

-  intensity  of  radiation  reflected  in  a 
direction  arc  cos  v  due  to  parallel  ravs 
of  incident  radiation  of  net  flux  -  in 
a  direction  arc  cos  u  (with  respect  to 
the  inward  directed  normal)  for  an 
isotropically  scattering  atmosphere  of 
optical  thickness  x  having  no  reflector 
at  the  bottom. 


The  functions  t(v,u,x)  and  T(v,u,x)  for  the  diffusely  trans¬ 
mitted  intensities  are  defined  similarly.  Also  introduce 

★  ★ 

the  functions  r  (v,u,x,A)  and  t  (v,u,x,A)  to  represent  the 
reflected  and  transmitted  intensities,  respectively,  for 
the  case  of  an  atmosphere  bounded  by  a  Lambert  law  re¬ 
flector  having  albedo  A. 

Reference  2  shows  that  the  functions  R  and  T  satisfy 
the  system  of  differential  equations: 


d 


R  ( v  ,  u ,  x ) 


R  { v  ,  u  ,  x ) 


F  v  ,  u  '  ,  x }  d u  '  ,-'u  1 


I  i  1 
|l*7i 

L  0 


u , x) dv 


0  %  v,  u  <  1 


X 
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R  ( v  ,  u  ,  0  )  -  o  , 


T  ( v  ,  u  ,  x ) 


T  ( v  ,  u  ,  x ) 


+  \  h  4-  4  /  R(v,u'  ,  x )  d u  '  'u' 


I  /  ,  L  I 

<  j  e  X'  U  +  ~  I  T  ( v  ’  ,u,x)dv*  /v*  j 


T  ( v  ,  u  ,  0  )  =  0 


The  parameter  ■  represents  the  air  edo  for  single  scattering. 

★  * 

To  obtain  equations  for  r  and  t  ,  let  r  (t>  ,x ,  A)  be  the 
constant  intensity  of  radiation  reflected  from  the  bettor 
surface  of  albedo  A,  the  incident  direction  being  arc  cos  u 
and  the  optical  thickness  being  x  [3].  Then  us  inn  the  con¬ 
servation  law  that  at  the  bottom  surface 


(upward  fluxi  A  •  (downward  flux)  , 


the  equation 


;  ~Ue 


( V  1  ,  U  ,  X  ) 


-  1  (  u  ,  X  ,  A  )  - 


•  i  n  /  ,  •  1  , ,  •  y  V 

4  i  I  ( u  ,  x  ,  A )  i  d 


!  o 


for  the  function  I  is  obtained.  Solving  for  I  vield. 
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I  (u  ,x ,  A) 


1  1 

1  -  A  /  /  R  (v  ’ 
0  0 


T  { v  '  ,  u ,  x )  d  v ' 


u ’ ,x) du ' dv ! 


(6) 


Viewing  the  function  I(u,x,A)  as  a  new  source  of  radi- 
tion  incident  on  the  bottom  of  the  atmosphere  leads  to  the 
relations 

•H 

t  (v,u,x,A)  =  T{v,u,x)/4v  (7) 


+ 


I  (u,x,M  ,  R  (v  ,u '  ,  x) 

TT  > 

0 


2ndu 1 


/ 


=  T (v,u,x) /4v 


+ 


1 

2 


I \ufx,A) 


1 

/  R(v,u' ,x)du'/v  , 
0 


(8) 


* 

r 


(v,u,x,A) 


=  R(v,u,x)/4v  +  I(u,x,A)e  X//v 


(9) 


+  I (u ,x  f A) 

IT 


I  T (v ,u 1 ,x) 
J  4v 
0 


2iTdu ' 


t 


-  R(v,u,x)/4v 


(10) 


+  I  (u,x,A) 


0 


T  (v , u " ,x)du '/v 
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The  above  equations  clearly  exhibit  the  fact  that 
knowledge  of  the  radiation  fields  for  the  nonplanetary 
problem  is  sufficient  for  obtaining  the  radiation  fields 
for  the  planetary  problem  with  a  Lambert-law  reflecting 
surface  [3]  . 
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III.  THE  COMPUTATIONAL  SCHEME 


For  the  purposes  of  numerical  calculation  of  the 

*  ★ 

quantises  r  (v,u,x,A),  t  (v,u,x,A),  and  the  reflected 
and  transmitted  fluxes,  finite  sums  replace  all  integrals 
appearing  in  Eqs .  (1),  (3),  (6),  (7),  and  (9).  If 
z^,Z2f...,ZjsJ  represent  the  nodes  and  ,w0 , . , . ,wN  the 
weights  of  an  N  point  quadrature  scheme,  the  above  dif¬ 
ferential  integral  equations  reduce  to  the  system  of 
ordinary  differential  equations 


Rij 


R,  .  (x)  =  - 


+  i- 

i  zi  2j 


R.  .  (x) 
1D 


N 


+  x  1  1  + 1  Jx  Rik(x)V'zk 


N 

l+5  l  Rn<x>wi/ 


k=l 


kj k'  k 


Ri j  (0)  =  0  , 


(11) 


(12) 


2-  T. . (x) 

dx  i] 


-  ~  Tij(x! 

l  J 


(13) 


+  X 


N 

1  + 1  Rik(x)Vzk 


N 

!"X/Zj  +  Kh  VX>Vzk 
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Ti.(0)  =  0  ,  i , j=l , 2 , . . . ,N  ,  (14) 

x  ^  0  . 


Here  the  convention 


j  (x)  =  Rfz^z^x)  ,  (15) 


T\_.  (x)  =  T  (zi  ,Zj  ,x)  , 


is  used.  Knowledae  of  the  functions  R. .  and  T, .  allows 

*  ID 

calculation  of  the  functions  r  and  t  via  the  formulas 


rij  (X(A) 


R, . (xJ/42^  +  I  ( z . ,x , A) 

J  ^  J 


(16) 


N 


X//zi  +  4  J  T..  (x)w,/z. 

2  lk  k'  l 


x>  0,  0  <  A  £  1  , 


t. . (x , A)  =  T. . (x)/4z. 

l]  '  i-}  1 


+ 


1 

2 


I  (Zj ,x , A) 


N 

J  R.k(x)wk/z1  , 


(17) 


x>0,0SASl. 


The  quantity  I(Zj,x,A)  is  calculated  by 
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I  ( z  j  ,  x ,  A )  = 


N 

z:°"X/Zi  +  5  (x)wk 


N  N 

1  -  A  .1  £, Rkm(x)“k“m 

k=l  m=l 


(18) 


The  reflected,  transmitted,  and  global  transmitted 
fluxes,  defined  by  the  equations 

1  * 

p(u,x,A)  =  2 tt  /  r  (z,u,x,A)z  dz  ,  (19) 

0 


1 

T  (u  ,  X  ,  A)  =  2 TT  I 


★ 

t  ( z ,u , x , A) z  dz  , 


0 


(2C) 


and 


Tg(u,x,A)  =  ?re  X/>U  +  t(u,x,A)  (21) 


are  computed  by  replacing  the  Integra"  '  with  sums  and  by 

★  * 

using  the  intensity  functions  r^  (x,A)  and  t^(x,A)  in  the 
formulas 


N 

P  ( z  .  ,  x ,  A)  =  2tt  l  r*.(x,A)z.w.  ,  (22) 

J  i=l 


N 

r  * 

x(z.,x,A)  =  2tt  l  t..(x,A)z.w.  , 
J  i=l  1  x 


S 


(23) 
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and 

xg(Zj,x,A)  =  it  6  X//Uj  +  T  (u  j  ,  x ,  A)  .  (24) 

The  following  steps  summarize  the  numerical  procedure: 

1)  Integrate  Eqs.  (11)  and  (13)  with  the  initial 
conditions  of  Eqs.  (12)  and  (14)  from  x  =  0  to 
x  =  x^,  the  desired  thickness. 

2)  At  x  =  x^,  for  a  fixed  value  of  A,  compute 
I(u,x,A)  from  Eq.  (18). 

3)  Using  the  computed  value  of  I(u,x,A)  and  the 
solutions  of  step  1,  calculate  the  reflected 
and  transmitted  intensities  from  Eqs.  (16) 
and  (17) . 

4)  Calculate  reflected  and  transmitted  fluxes 
from  Eqs.  (21)  and  (22)  using  previously 
computed  intensities  of  step  3. 

Note:  Since  the  parameter  A  occurs  only  in  the  calculation 

of  I(u,x,A),  the  set  of  intensities  and  fluxes  for  a  range 
of  A  values  may  be  obtained  by  recalculating  only  I(u,x,A) 
for  the  desired  values  of  A. 


-10- 


IV.  NUMERICAL  RESULTS 


The  basic  numerical  calculation  consisted  of  producing 
the  reflected  and  transmitted  intensities  and  fluxes  for 
surface  albedos  A  =  0 , .  1 , .  2  , . .  .  ,  1 . 0  ,  optical  thickness  0-100, 
and  conservative,  isotropic  scattering  (A  =  1),  The  integra¬ 
tion  step  size  used  was  A  =  .005  with  Gaussian  quadrature  of 
order  N  =  7.  A  check  calculation  involving  changing  the  step 
size  to  A  =  .01  and  varying  che  order  of  the  quadrature  from 
N  =  3  to  N  =  5  was  run  to  a  thickness  of  1  resulting  in,  at 
most,  a  change  of  one  unit  on  the  fourth  significant  figure. 
All  calculations  were  performed  on  a  CDC  6600  computer  using 
a  fourth-order  Adams-Moulton  predictor-corrector  integration 
scheme.  Execution  time  for  the  basic  calculations  was  about 
7  min. 

One  primary  objective  of  the  computational  study  was 
to  evaluate  the  differences,  if  any,  in  the  reflected  and 
transmitted  intensities  and  fluxes  when  both  Rayleigh  and 
isotropic  scattering  laws  were  assumed.  Reference  1  has 
examined  the  case  of  Rayleigh  scattering  by  solving  the 
problem  via  singular  integral  equations.  Comparison  of 
Figs.  1,  2,  and  3  (depicting  reflected  diffusely  trans¬ 
mitted,  and  global  transmitted  fluxes)  with  the  analogous 
graphs  in  Kahle's  paper  leads  to  this  conclusion:  as  far 
as  reflected  and  transmitted  fluxes  are  concerned,  virtually 
no  quantitative  difference  exists  between  the  reflected  and 
transmitted  fluxes  tor  the  two  scattering  laws  for  any 
optical  thickness  in  the  range  0-100  and  any  surface  albedo 
0  £  A  £  1  for  normal  incidence. 

Tables  1  to  9  present  the  intensities  and  fluxes  of  the 
reflected  and  diffusely  transmitted  radiation  for  slab  thick¬ 
nesses  .2,  10,  and  100  and  surface  albedos  A  =  0,  .5,  and 
1.0.  These  tables  have  been  excerpted  from  the  main  calcu¬ 
lations  *-hat  output  intensities  and  fluxes  at  61  optical 


thicknesses  in  the  range  0-100  for  A  =  0 . 0 , 0. 1, 0 . 2 , . . .  ,  1. 0 . 
The  incident  angle  is  constant  across  a  row  in  the  tables. 
It  takes  on  one  of  eight  different  angles  as  indicated. 

The  tables  give  intensities  for  seven  outgoing  angles,  and 
list  fluxes  in  the  last  column. 


REFLECTED  AND  DIFFUSELY  TRANSMITTED  INTENSITIES  AND  FLUXES  FUR 
SLAB  WITH  CONSERVATIVE  ISOTROPIC  SCATTERING,  THICKNESS  0.?, 


Table  2 

REFLECTED  AND  DIFFUSELY  TRANSMITTED  INTENSITIES  AND  FLUXES  FOR 
SLAB  WITH  CONSERVATIVE  ISOTROPIC  SCATTERING,  THICKNESS  0.2, 
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V.  COMPUTATIONAL  CHECKS 


Checks  on  the  numerical  results  may  be  divided  into 
two  categories:  1)  internal  checks  on  the  consistency  of 
the  numerical  scheme;  and  2)  external  checks  using  other 
sources  and  independent  methods. 

ontioned  previously,  the  internal  checks  consist 
of  changing  the  integration  step  size  from  A  =  .005  to 
A  =  .01  and  varying  the  order  of  the  quadrature  scheme  from 
N  =  7  to  N  =  3  and  5.  In  all  cases,  the  optical  thickness 
ranged  from  0  to  1;  and  the  results  changed  by  no  more 
than  one  unit  in  the  fourth  significant  figure. 

The  first  external  check  considered  is  restricted  to 
the  case  A  =  0.  Reference  2  presents  tables  similar  to 
Tables  1  through  9  for  optical  thicknesses  1-50.  In  all 
cases,  the  current  results  agree  perfectly  with  those 
earlier  calculations. 

A  second  check,  for  the  case  A  =  1,  normal  incidence, 
is  the  conservation  relation  that  the  reflected  flux  must 
equal  the  incident  flux.  At  normal  incidence,  the  re¬ 
flected  flux  is  3.1416  and  at  an  angle  of  incidence  of  60° 
it  is  1.5708  -  n / 2 .  As  Tables  3,  6,  and  9  exhibit,  this  is 
the  case  for  all  optical  thicknesses. 

An  independent  check  is  the  calculation  of  reflected 
and  transmitted  intensities  by  generalized  X  and  Y  func¬ 
tions  as  outlined  in  Ref.  3.  Both  the  cases  u  =  0.0, 
v  =  cos  88.5°,  and  u  =  0.5,  v  =  cos  88.5°  for  A  =  .8  and 
optical  thickness  x  =  0.8  have  been  computed  by  this  method 
giving  results  that  agree  exactly  with  those  calculated  by 
the  initial-value  procedure. 

The  final  check  compares  qualitatively  the  results  of 
the  initial-value  method  to  Kahle's  results  [1).  Even  though 
Kahle  assumed  Rayleigh  scattering,  both  sets  of  calculations 
should  be  qualitatively  the  same.  Comparing  Figs.  1  through 
3  to  the  analogous  figures  in  Kahle's  paper  bears  out  this 
conjecture . 
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An  interesting  property  of  the  solution  of  the  planetary 
problem  with  a  Lambert- law  reflector  is  obtained  by  comparing 
this  solution  with  that  presented  in  Ref.  4  for  the  case 
when  the  reflecting  surface  is  a  perfect  specular  reflector. 
It  is  observed  for  the  conservative  case  A  =  1,  A  =  1  that 
the  transmitted  and  reflected  fluxes  are  virtually  the  same 
(to  three  figures)  for  all  optical  thicknesses  ^  3.  Thus, 
even  the  properties  of  the  reflecting  surface  have  little 
effect  on  fluxes  for  moderately  thick  media. 
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